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Abstract 

Dirac Hamiltonian is scaled in the atomic units h = m = 1, which allows us to take the non- 
relativistic limit by setting the Compton wavelength A — >■ 0. The evolutions of the spin and 
pseudospin symmetries towards the non-relativistic limit are investigated by solving the Dirac 
equation with the parameter A. With A transformation from the original Compton wavelength 
to 0, the spin splittings decrease monotonously in all spin doublets, and the pseudospin splittings 
increase in several pseudospin doublets, no change, or even reduce in several other pseudospin 
doublets. The various energy splitting behaviors of both the spin and pseudospin doublets with A 
are well explained by the perturbation calculations of Dirac Hamiltonian in the present units. It 
indicates that the origin of spin symmetry is entirely due to the relativistic effect, while the origin 
of pseudospin symmetry cannot be uniquely attributed to the relativistic effect. 
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It is well known that the spin and pseudospin symmetries play critical role in the shell 
structure and its evolution. The introduction of spin-orbit potential made the single- 
particle shell model can well explain the experimentally observed existence of magic num- 
bers for nuclei close to the valley of ^-stability To understand the near degener- 
acy observed in heavy nuclei between two single-particle states with the quantum numbers 
(n — 1, 1 + 2, j = I + 3/2) and (n, l,j = l + 1/2), the pseudospin symmetry (PSS) was intro- 
duced by defining the pseudospin doublets {n = n — 1,1 = I + l,j = I ±1/2) 3, |4[, which 
has explained numerous phenomena in nuclear structure including deformation 5|, superde- 
formation Jfj], identical bands [7], and magnetic moment jgj]. Because of these successes, 
there have been comprehensive efforts to understand their origins as well as the breaking 
mechanisms. For the spin symmetry (SS), the spin-orbit potential can be obtained natu- 
rally from the solutions of Dirac equation. Thus, the SS can be regarded as a relativistic 
symmetry. For the PSS, its origin has not been fully clarified until now. It is worth re- 
viewingsome of the major progresses in understanding the underlying mechanism of PSS. 
In Ref. [91], a helicity unitary transformation of a non- relativistic single-particle Hamiltonian 
was introduced to discuss the PSS in the non-relativist ic harmonic oscillator. The particular 
condition between the coefficients of spin-orbit and orbit-orbit terms was indicated in the 
corresponding non-relativistic single particle Hamiltonian for the requirement of PSS. The 
same kind of unitary transformation was considered in Ref. 10(, where the application of the 
helicity operator to the non-relativistic single-particle wavefunction maps the normal state 
(l,s) onto the pseudo-state (7, s), while keeping all other global symmetries. A substantial 
progress was achieved in Ref. rjj, where the relativistic feature of PSS was recognized. The 
pseudo-orbital angular momentum I is nothing but the orbital angular momentum of the 
lower component of Dirac spinor, and the equality in magnitude but difference in sign of the 
scalar potential S and vector potential V was suggested as the exact PSS limit. Meng et al. 
showed that exact PSS occurs in the Dirac equation when the sum of the scalar S and vector 
V potentials is equal to a constant [12j |. Unfortunately, the exact PSS cannot be met in real 



nuclei, much effort has been devoted to the cause of splitting. In Ref. |13l-ll5|. it was pointed 
out that the observed pseudospin splitting arises from a cancellation of the several energy 
components, and the PSS in nuclei has a dynamical character. A similar conclusion was 
, In addition, it was noted that, unlike the spin symmetry, the pseu- 
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dospin breaking cannot be treated as a perturbation of pseudospin-symmetric Hamiltonian 
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18f | . The non-perturbation nature of PSS has also been indicated in Ref.[19|. Regardless of 
these pioneering studies, the origins of the spin and pseudospin symmetries have not been 
fully understood in the relativistic framework. Recently, we have checked the PSS by use 
arity renormalization group and shown explicitly the relativistic origin of this 



of the simi 

symmetry [20|. However, the dependence of the quality of PSS on the relativistic effect has 
not been checked until now. In this paper, we study the evolution of the spin and pseudospin 
symmetries from the relativistic to the non-relativistic to explore the relativistic relevance 
of this symmetries. 

Dirac equation of a particle of mass m in external scalar S and vector V potentials is 
given by 

H = ca ■ p + f3 (mc 2 + S) + V, (1) 

where a and /3 are the usual Dirac matrices. For a spherical system, the Dirac spinor ip has 
the form 

^ = 1 | iG nK (r) (f) Kmj (0, if) \ ^ 
r Y F nK (r) a ■ r(/) Kmj (i?, </?) J 
where n is the radial quantum number, and rrij is the projection of angular momentum on 
the third axis, k = ±(j + 1/2) with — for aligned spin (si/ 2 , P3/2, etc.), and + for unaligned 
spin (pi/2, d 3 /2, etc.). Splitting off the angular part and leaving the radial functions satisfy 
the following equation 



mc 2 + S(r) -cf r + f \(G(r)\ (G(r) 
c£ + f -mc 2 + A(r) ) [ F (r) I ~~ ' [ F (r) 



(3) 



where £(V) = V(r) + S(r) and A(r) = V(r) — S(r). Based on Eq.()3]), a lot of work has been 



carried out to check the origins of the spin and pseudospin symmetries [Hi, I13H15L Il8l . Il9 |. 
Although they are recognized as the symmetries of Dirac Hamiltonian, it is still not very 
clear the important role of relativistic effect. In order to explore the relativistic effects of this 
symmetries, the atomic units H — m — 1 are adopted instead of the conventional relativistic 
units H = c = 1 in the present system. For simplicity, the operator H is measured in unit of 
the rest mass, mc 2 . Then the equation ([3]) is presented as 



1 + A 2 S A (-£ + *) 
Hi + 7) - 1 + A2A 
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where the Compton wavelength A = h/mc — 1/c. In this units, the result in the non- 
relativistic limit can be obtained in a very simple, intuitive, and straightforward manner by 
taking the speed of light c — > oo or the Compton wavelength A — > 0, which is not possible 
in the latter units since c = 1. 

In order to investigate the evolution from the relativistic to the non-relativistic, A is 
regarded as a parameter and the original Compton wavelength A = h/mc is labelled as 
Ao- The relativistic result corresponds to the solution of Eq.fll]) with A = Ao- The result 
in the non-relativistic limit can be obtained from Eq.(j3]) by setting A — > 0. Thus, the 
evolution from the relativistic to the non-relativistic can be checked by transforming A 
from Ao to 0. Then, the relativistic effects of the spin and pseudospin symmetries can be 
investigated by extracting the energy splittings between the spin or pseudospin doublets, 
and this symmetries develop toward the non-relativistic limit can be checked, and vice versa. 

In order to make this clear, we have solved Eq.fll]) for a Woods-Saxon type potential for 
S(r) and A(r), i.e., S(r) = E /(as,Ts,^) and A(r) = A /(aA, ?"a, r) with 

/(a ,r ,r) = r-. (5) 

1 + ex P (^r) 

The corresponding parameters are determined by fitting the energy spectra from the RMF 
calculations for 208 Pb (to see Ref. 2lj). The energy spectra of Eq.(jl]) are calculated by 
expansion in harmonic oscillator basis. 

The single particle energy varying with the parameter A is displayed in Fig.l, where it 
can be seen that the energy decreases monotonously with A decreasing for all the levels 
available. The trend of energy with A is towards the direction of non-relativistic limit. 
With the decreasing of A, the calculation is closer to the non-relativistic result. When A is 
reduced to A/Ao = 0.1, the solution of Eq.(jl]) is almost same as the non-relativistic result. 
Furthermore, for the different single-particle states, the sensitivity of energy to A is different. 
For the spin unaligned states, the decreasing of energy is faster than that for the spin aligned 
states, which leads to the energy splittings of the spin doublets reduce with A decreasing. 
When A is reduced to A/Ao = 0.1, the spin-orbit splittings almost disappear for all the spin 
doublets. These indicate that the spin symmetry becomes better as A decreases, and the 
spin symmetry breaking is entirely due to the relativistic effect. 

To better understand the preceding claim, the energies in several A values are listed in 
Table I for all single particle levels. For comparison, Table I does also display the data 
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of the non-relativistic calculations (the last column), which are obtained by solving the 
Schrodinger equation Hip(r) = Eip(r) with H = — ^ — ^I^M + S(r). From Table I, 
it can be seen that the relativistic spin-orbit splitting (A = Ao) is considerably large. This 
splitting decreases with the decreasing of A. When A/Ao = 0.001, the energy of the spin 
unaligned state in conjunction with that of the spin aligned state degenerates to the non- 
relativistic result. These indicate that Eq.fll]) reproduces well the process of development 
from the relativistic to the non-relativistic, and both results of the relativistic and non- 
relativistic can be obtained well from Eq.(j3]) with an appropriate value of A. Hence, the 
relativistic effects of the spin and pseudospin symmetries can be checked from the solutions 
of Dirac equation with the parameter A. 

In order to recognize clearly the relativistic effect of spin symmetry, the energy splittings 
of spin doublets varying with A are plotted in Fig.2, where it is shown that the energy 
splittings decrease monotonously with A reducing for all the spin partners. When A is 
reduced to A/Ao = 0.1, the energy splittings of all the spin doublets are almost reduced to 
zero. The detailed observation shows that the energy splittings are more sensitive to A for 
the states with higher orbital angular momentum in the same radial quantum number. For 
the states with the same orbital angular momentum, the energy splittings appear crosses in 
the different radial quantum number. These reflect that the relativistic sensitivity is different 
for the states with different quantum numbers. When A is reduced to zero, the spin-orbit 
splittings disappear for all the spin partners, the non-relativistic results are obtained in 
excellent agreement with those from the solutions of Schrodinger equation. Namely, the 
spin-orbit splitting arises completely from the relativistic effect ,and can be treated as a 
perturbation of spin-symmetric Hamiltonian as indicated in Ref. 18]. 

Different from the spin symmetry, the relativistic origin of pseudospin symmetry is more 
complicated. In Fig. 3, we display the energy splittings of pseudospin doublets varying 
with the parameter A. From there, it can be observed that the energy splittings increase 
significantly with A decreasing for the pseudospin partners (2g 9 / 2 ,lin/2), (2/7/2,1^9/2), and 
(3^3/2,2/5/2). Especially for (2g Q / 2 ,lin/2), the increasing of energy splitting is very obvious. 
For the doublets (2^5/2,1(77/2), the increasing of energy splitting with A decreasing is relatively 
small. When A/A decreases below than 0.6, the energy splitting goes toward a stable value. 
The same phenomenon also appears in the doublet (3si/ 2 ,2<i3/2). However for the pseudospin 
doublets (2^3/2,1/5/2) and (2^1/2,1^3/2), an opposite evolution of energy splitting with A is 
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disclosed. It shows the origin of pseudospin symmetry is more complicated than that of spin 
symmetry. The pseudospin splitting cannot be attributed uniquely to the relativistic effect. 
The quantum number of single-particle states and the shape of potential make important 
influence on this symmetry. 

In order to better understand the relativistic effects of the spin and pseudospin sym- 
metries, we expand perturbatively the Dirac Hamiltonian in Eq.Q to analyze the effects 
of each higher order term on the energy splitting behaviors of both spin and pseudospin 
doublets. Following Ref. 20|], for Dirac particle, the expanded Hamiltonian up to the order 
1/m 3 is 

ff = E(r) + ^-A 2 l( V - S'4-) - A 2 ^ + A 4 A ( V - 2S'4-) + \*«? A ' 



2m 2m? V dr J r Am 2 2m? V dr I r 2m? 



2 E" 2 p 4 4 S" -2£'A' + 4S£" 



+A 2 -^ - A r - A (6) 

8m 2 8m 6 16m 6 



where p 2 = —4^ + K t ■ Based on the same considerations as Ref. 



H is 



dr 2 

2 - " 1 ' " " "■ x - 2 k A' 

r 4m 2 ' 

' 2 _OV A' 



decomposed into the eight components: £ (r) + 2—, — A 2 ^- (Sp 2 — S'-^p), — A 2 ^ 
+\^(Sp 2 -2S>i), + A*f0, +A 2 ^, -A 2 ^- , - >' a -*£^'' , which are respec- 
tively labelled as Oi, O2, • ■ ■ ,0$. 0\ corresponds to the Hamiltonian in the non- relativistic 
limit, i.e., the Schrodinger part of H . 02(04) is the dynamical term relating to the order 
l/m 2 (l/m 3 ). 03(05) is the spin-orbit coupling corresponding to the order l/m 2 (l/m 3 ). The 
eigenvalues of H are calculated with the fully same E(r) and A(r) as that in calculating the 
exact solutions of Eq.fll]). 

For recognizing the relativistic effect of SS, we analyze the reason why the energies of the 
spin unaligned states decrease faster than those of the spin aligned states when A decreases. 
As an illustrated example, we display the energy splittings of every component Oj(i = 
1,2, ••• ,5) varying with A for the spin doublets (lpi/2, IP3/2) and (1*77/2, 1^9/2) in FigJU 
where we neglect the results of 6 , Oy,and 8 because their contributions to the energy 
splitting are minor and do not influence on the total energy splitting behavior with A. From 
FigJU h can be seen that the contributions of all the Oj(i = 2, 3, 4, 5) to the energy splittings 
between the spin unaligned states and the spin aligned states are positive, and the positive 
energy splittings decrease with A decreasing. It is for this reason that the energies of the 
spin unaligned states decrease faster than those of the spin aligned states with A decreasing. 
Compared with O3 (the spin-orbit coupling corresponding to the order 1/m 2 ), and the 
contributions of 2 , 4 , and 5 to the spin energy splittings are relatively minor. The total 
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energy splittings are dominated by the contribution of O3 when A is sufficiently small. This 
means that, as the relativistic effect becomes weak, the spin splittings are almost entirely 
due to the spin-orbit coupling. For the different spin partners, the energy splitting behaviors 
with A are same except for the extent of splittings, as displayed in Fig JH for the spin doublets 
(lpi/2) I-P3/2) an d (1^7/2, 1^9/2)- These indicate that the spin symmetry origins completely 
from the relativistic effect, and possesses the perturbation attribute claimed in Ref.[18]. To 
understand the relativistic effect of PSS, we analyze the cause of the various energy splitting 
behaviors of pseudospin doublets. In Figj5j we show the energy splittings of each component 
Oi(i = 1, 2, • ■ ■ ,5) varying with A for the pseudospin partners (2si/ 2 , 1^3/2) and (2/7/2, 1/2-9/2)- 
From there, it can be seen that the pseudospin energy splittings caused by the Schrodinger 
part of H are dominated. This splittings are reduced by the contribution of spin-orbit 
coupling, and added by the contribution of dynamical terms. For the pseudospin partner 
(2si/2, 163/2), with the decreasing of A, the contribution of the pseudospin breaking (the 
dynamical terms) declines faster than that of the pseudospin improvement (the spin-orbit 
coupling), which results in better PSS when A decreases. However for the (2/7/2, I/19/2), the 
energy splittings caused by the pseudospin breaking varying with A are relatively slower than 
that by the spin-orbit coupling, which leads to the PSS becomes worse with A decreasing. 
These cause the different energy splitting behaviors of pseudospin doublets with A. Hence, 
the pseudospin splitting can not be regarded as a perturbation in agreement with the claim 
in Ref.[18|. 

In addition to the energy splittings associated with the relativistic effects, the wave 
function splittings between the (pseudo)spin doublets are also associated with the relativistic 
effects. An illustrated example is displayed in Fig|6j where the upper component of Dirac 
spinor for the states 1*77/2,9/2 is depicted in several A values. From FigJHl it can be seen that 
the wavefunction splitting of spin doublet is obvious for a relativistic particle (A/Ao = 1). 
With the development towards the non-relativistic direction (to reduce A), the wavefunction 
splitting of spin doublet decreases, which is in agreement with the case of level splitting. For 
the pseudospin symmetry, the lower component of Dirac spinor for the pseudospin doublet 
(2*79/2, H11/2) is drawn in FigJTJin several A values. The wavefunction splitting of pseudospin 
doublet is obvious when A/Ao = 1. Different from the spin splitting, we can not see that 
the pseudospin splitting reduces with A decreasing, which is consistent with the case of level 
splitting. 
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In summary, Dirac Hamiltonian is scaled in the atomic units h = m = 1, which allows 
us to take the non-relativist ic limit by setting the speed of light c — > oo or the Compton 
wavelength A — » 0. The evolution towards the non-relativistic limit is investigated from 
the solutions of Dirac equation by a continuous transformation of the parameter A. The 
solutions of Dirac equation corresponding to A = h/mc and A = represent respectively the 
relativistic result and that in non-relativistic limit. To transform the parameter A from h/mc 
to 0, the solutions of Dirac equation show the evolution from the relativistic to the non- 
relativistic limit. The relativistic effects of the spin and pseudospin symmetries are checked 
from the solutions of Dirac equation with the parameter A. It shows the spin splittings 
decrease monotonously with A reducing for all the spin partners. When A is reduced to 
zero, the spin-orbit splittings disappear, which is in agreement with the result in the non- 
relativistic calculations. For the pseudospin symmetry, the energy splittings increase in 
several partners, no change, or even decrease in another some partners. Compared with 
the spin symmetry, the origin of pseudospin symmetry is more complicated, and cannot be 
attributed uniquely to the relativistic effect. The quantum number of single-particle states 
and the shape of potential make important influence on this symmetry. By the perturbation 
calculations of Dirac Hamiltonian, the various energy splitting behaviors of both spin and 
pseudospin doublets with A are explained, which origins from the different contributions of 

n 

each component to the energy splittings. The result supports the claim in Ref . [18j , the spin 
splitting can be treated as a perturbation, while the pseudospin splitting can not be regarded 
as a perturbation. The same conclusion can also be obtained from the wavefunction. 
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KJ2009A129, and the 211 Project of Anhui University. 



[1] O. Haxel, J.H.D. Jensen, and H.E. Suess, Phys. Rev. 75, 1766 (1949). 

[2] M. Goeppert-Mayer, Phys. Rev. 75, 1969 (1949). 

[3] K.T. Hecht and A. Adler, Nucl. Phys. A 137, 129 (1969). 



8 (March 3, 2013) 
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dynaml, spin-orbl, dynam2, spin-orb2, and the total energy splitting is labelled as total. 
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FIG. 5: (Color online) The same as FigJH but for the pseudospin energy splittings. 
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